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1 Introduction and main results 

Triply orthogonal coordinate systems have attracted the attention of math- 
ematicians and physicists for almost two hundred years now. Particular 
examples of them were already used by Leibniz and Euler to evaluate mul- 
tiple integrals in canonical coordinates, and later Lame and Jacobi carried 
out calculations in analytical mechanics with the help of the famous elliptic 
coordinates. The first general results on the geometry of triply orthogo- 
nal systems date back to the 19th century - like the famous theorem of 
Dupin, saying that coordinate surfaces intersect along curvature lines. The 
Lame equations which analytically describe the triply orthogonal systems 
were studied in detail by Bianchi Bi2j and Darboux [Da . 



Recently orthogonal systems came back into the focus of interest in 
mathematical physics as an example of an integrable system. Zakharov 
[Z] has shown how the Lame equations can be solved by the <9-method 
and constructed a variety of explicit solutions with the help of the dress- 
ing method. Algebro-geometric solutions of the Lame equations were con- 
structed by Krichever jKj. The recent interest to the orthogonal coordinate 
systems is in particular motivated by their applications to the theory of the 
associativity equations [Du . 
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and the DFG Research Center "Mathematics for Key Technologies" (FZT 86) in Berlin. 
^Supported by the SFB 288 "Differential Geometry and Quantum Physics". 
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Figure 1: An elementary hexahedron of a discrete orthogonal system 



The question of proper discretization of the classical models in differential 
geometry became recently a subject of intensive study (see, in particular, 
|BP| IDS2|. IKSj ) . Indeed one can suggest various models in discrete geometry 
which have the same continuous limit and nevertheless have quite different 
properties. For a great variety of geometric problems described by integrable 
equations it was found that integrable discretizations (i.e. the discretizations 
preserving the integrability of the underlying nonlinear system) preserve the 
characteristic geometric properties of the problem. Moreover, it turns out 
that this discretization can be described in pure geometric terms. 

Such a discretization of the triply-orthogonal coordinate systems was first 
suggested by one of the authors |Boj . It is based on the classical Dupin the- 
orem. Since the circular nets are known |MPS[ IN] to correspond to the 
curvature line parametrized surfaces, it is natural to define discrete triply or- 
thogonal systems as maps from 1? (or a subset thereof) to M. 3 with all elemen- 
tary hexahedrons lying on spheres. The neighboring spheres intersect along 
circles which build the coordinate discrete curvature line nets (see FigOJ). 
Doliwa and Santini [CDS j made the next crucial step in the development 
of the theory. They considered discrete orthogonal systems as a reduction 
of discrete conjugated systems |DSlj . generalized them to arbitrary dimen- 
sion and proved their geometric integrability based on the classical Miguel 
theorem |Bej . The latter claims that provided the seven (black) points lie 
on circles as shown in Fig^ the three dashed circles intersect in a com- 
mon (white) point. This implies that a discrete triply-orthogonal system is 
uniquely determined by its three coordinate circular nets (as shown in FigEI) • 
This is a well-posed initial value problem for triply orthogonal nets. 

Later on discrete orthogonal systems were treated by analytic methods of 
the theory of solitons. Algebro-geometric solutions |AKVj as well as the d- 
method |DMSj were discretized preserving all the symmetries of the smooth 
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Figure 2: Initial data for a smooth Figure 3: Initial data for a discrete 
triply orthogonal system. triply orthogonal system. 



case. 

There is a common belief that the smooth theory can be obtained as a 
limit of the corresponding discrete one, which can be treated also in this 
case as a practical geometric numerical scheme for computation of smooth 
surfaces. Moreover, this discrete master theory has an important advantage: 
surfaces and their transformations appear exactly in the same way as sublat- 
tices of multidimensional lattices. This description immediately provides 
the statements about the permutability of the corresponding (Backlund- 
Darboux) transformations, - an important property which is nontrivial to 
see and sometimes even difficult to check in the smooth theory. 

Until recently there were no rigorous mathematical statements supporting 
the observation about the classical limit of the discrete theories. The first step 
in closing this gap was made in our paper [BMSJ where general convergence 
results were proven for a geometric numerical scheme for a class of nonlinear 
hyperbolic equations. The first geometric example covered by this theory 
was smooth and discrete surfaces with constant negative Gaussian curvature 
described by the sine-Gordon equation. 

In the present paper we formulate initial value problems and prove the cor- 
responding convergence results for conjugate and orthogonal coordinate sys- 
tems, which are described by three-dimensional partial differential/ difference 
equations. One should note here that despite the known geometric dis- 
cretization, a convenient analytic description of discrete orthogonal systems 
converging to the classical description of the smooth case was missing. In 
particular this was the problem with discrete Lame equations (cf. [AKVj). 
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The Clifford algebra description of smooth and discrete orthogonal systems 
suggested in |BH| turned out to be optimal for our purposes. This descrip- 
tion provides us with the proper discrete analogs of the continuous quantities 
and their equations. 

In particular we introduce the discrete Lame system which is the start- 
ing point for further analytical investigations. We write these equations in 
the form of a discrete hyperbolic system and formulate the corresponding 
Cauchy problem. As a by-product, we rewrite the classical Lame system 
in the hyperbolic form and pose a Cauchy problem for classical orthogonal 
coordinate systems, which seems to be new. 

The central theorems of this paper are on the convergence of discrete 
orthogonal coordinate systems to a continuous one. In particular we prove 
the following statement. 

Theorem 1 Assume three umbilic free immersed surfaces ?i ; 3 ^ and 9^ in 
R 3 intersect pairwise along their curvature lines, Xi = 9^ n 33 etc. Then: 

1. On a sufficiently small box B = [0,r] 3 C R 3 , there is an orthogonal co- 
ordinate system x : B — > R 3 ; which has 9i, 9 r 2 and 9^ locally as images 
of the coordinate planes, with each curve Xi being parametrized locally 
over the i-th axis by arc-length. The orthogonal coordinate system x is 
uniquely determined by these properties. 

2. Denote by B e = B H (eZ) 3 the grid of mesh size e inside B. One 
can define a family {x e : B e — > 1R 3 } of discrete orthogonal systems, 
parametrized by e > 0, such that the approximation error decays as 

sup - < Ce. 

The convergence is C°° : all partial difference quotients of x e converge 
with the same rate to the respective partial derivatives of x. 

We also answer the question of how to construct the approximating fam- 
ily {x e }: One solves the discrete Lame system with discrete Cauchy data 
determined by 3i, 9^ and 9^. The solution x e is then calculated in a process 
that consists of 0(e _3 )-many evaluations of the expressions in the discrete 
Lame equations. The scheme can be - and has been - easily implemented 
on a computer, it needs a precision of e 2 for the calculations and is robust 
against C°°-deviations of order e in the initial data. 

Actually, the results presented here go further: 

• We do not restrict ourselves to the case of three dimensions, but allow 
arbitrary dimensions, possibly different for domain and target space. 
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In particular, we prove the corresponding approximation results for 
curvature line parametrized surfaces. 

• One can keep some of the directions of the orthogonal system discrete 
while the others become continuous. This way, the classical Ribaucour 
transformation is obtained very naturally. The approximation result 
holds simultaneously for the system and its Ribaucour transforms. 

• The results about the permutability of the Ribaucour transformations 
are elementary in the discrete setup. Our approximation theorem im- 
plies the corresponding claim for classical orthogonal systems (Theorem 
IHJ). The first part of this theorem is due to Bianchi |Bilj . the second 
part was proved in |GTj using Bianchi's analytic description of the 
Ribaucour transformations. 

• In addition to orthogonal systems, we also treat, in a completely anal- 
ogous manner, conjugate nets. We pose a Cauchy problem for them 
and show that any continuous conjugate net, along with its transfor- 
mations, can be C°°-approximated by a sequence of discrete conjugate 
nets. The latter are quadrilateral lattices introduced by Doliwa and 
Santini [DS1. 

The paper is organized as follows. We start in Sect. El with a short pre- 
sentation of the concepts developed in our recent paper [BMS j5 stating the 
main result about an approximation of solutions of a hyperbolic system of 
partial differential equations by solutions of the corresponding partial dif- 
ference system. In Sect. |3] we apply this theorem to conjugate systems. 
We formulate a Cauchy problem for smooth and discrete conjugate systems, 
define the Jonas transformation and prove the convergence. Sect. [His dedi- 
cated to the description of orthogonal systems in Mobius geometry. Sect. El 
contains the approximation results for orthogonal systems, a part of which 
was summarized above. 

2 The approximation theorem 

We start with a study of hyperbolic systems of partial difference and partial 
differential equations. A general convergence theorem is formulated, stating 
that solutions of difference equations approximate, under certain conditions, 
solutions of differential equations, provided the approximation of equations 
themselves takes place. 
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2.1 Notations 



For a vector e = (ei, . . . , ejvf) of M non-negative numbers, we define the 
lattice 

T = (eiZ) x ... x (e M Z) c K M (1) 

We refer to as the mesh size of the grid IP 6 in the z-th direction; a direction 
with a vanishing mesh size = is understood as a continuous one, i.e. in 
this case the corresponding factor e{L in (0) is replaced by R. For instance, 
the choice e = (0, . . . , 0) yields T = R M . We use also the following notation 
for n-dimensional sublattices of D 56 (n < M): 

n , = {e = ■ • • , 6*) e ? e l & = o if < Oi, . . . , *n» . (2) 

Considering convergence problems, we will mainly deal with the following 
two situations: 

1. €i = e for all i — 1, . . . , M (all directions become continuous in the limit 

2. €i = e for i = 1, . . . ,m; €i = 1 for i = m + 1, . . . , M (the last M — m 
directions are kept discrete in the limit). 

Thus, we always have one small parameter e > only. In our geometric prob- 
lems, the first situation corresponds to the approximation of conjugate and 
orthogonal systems themselves, while the second situation corresponds to the 
approximation of such systems along with their transformations. The under- 
lying difference equations for discretizations of conjugate (resp. orthogonal) 
systems, and for their transformations will be the same, the only distinction 
being in the way the continuous limit is performed. We will distinguish be- 
tween these two situations by setting M = m + m', where m! — and in the 
first situation, while w! > in the second one. 

Our considerations are local, therefore the encountered functions are usu- 
ally defined on compact domains, like the cube 

B(r) = {£ G R M | < £ < r for i < M] (3) 

in the continuous context, and the following compact subsets of IP 6 in the 
lattice context. In the situation 1 (m' = 0): 

B e (r) = {£ G ? e | < £ < r for i < M} ; (4) 

in the situation 2 [m! > 0): 

B e (r) = {£ G T | < ii < r for 1 < i < m, & G {0, 1} for m < % < M} . 

(5) 
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We will often use the notation 



B £ (r) = {£ G (eZ) m | < & < r for 1 < i < m} , (6) 

so that in the latter case B e {r) = Bq(t) x {0, l} m '. 

For lattice functions u : T e i— > X with values in some Banach space X, we 
use the shift operators Tj and the partial difference operators <5j = (Tj — 1)/^ 
defined as 

(r i u)(0 = «(^ + e i e i ), = ~ + ~ 5 ( 7 ) 

here is the ith unit vector. For = 0, define = 1 and 5i = di. Multiple 
shifts and higher order partial differences are written with the help of multi- 
indices a = (ai, . . . , ckm) £ N M : 

7 .a =r ai... 7 .« M> S a = 5^---S a M M , \a\=a 1 + ... + a M . (8) 
We use the following C^-norms for functions u : B e (r) — > X: 

|M|, = sup{|^0|x|£e^(r-eM); M < I (9) 

if m' > 0, only the multi-indices a with at = (« > m) are taken into 
account in the last formula. Any function defined on B t (r) with e > has 
finite C l - norms of arbitrary order l\ on the other hand, if u : B(r) i— > X 
is a function of a continuous argument which does not belong to the class 
C (B(r)), then for its restriction to the lattice points, u e — u\B<^(r), the norm 
will, in general, diverge for e — > 0. 

Definition 1 VFe sa?/ £/ia£ a family of maps {u e : B e (r) — > X}o< e <ei 0(e)- 
convergent in C°°(B(r)) to a smooth limit u : -B(r) — > X ; if there are finite 
real constants such that 

\\u e - u\ B *{r) \\e < K e e. 
for any integer £ > and all positive e < r/£. 



2.2 Approximation for hyperbolic equations 

Being interested in the convergence of solutions of difference equations to- 
wards the solutions of differential ones, we restrict our attention to two ver- 
sions of the limit behaviour of the lattice J 6 described above, i.e. m' — (all 

6j = e) or m' > (e* = e for 1 < i < m, and ej = lform + l<-i<M = 
m + ml). 
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Let X — Xi x ■ ■ ■ x X K be a direct product of K Banach spaces, equipped 
with the norm 

KOIx= max K(£)|x fc . (10) 

l<k<K 

Consider a hyperbolic system of first-order partial difference (e > 0) or partial 
differential (e = 0) equations for functions u : D 56 i— > X: 

Vfe = /fc», 1<A:<^, je£(fc). (ll) 

Thus, for each = 1, there is defined a subset of indices E(k) C 

{1, . . . , M}, called the set of evolution directions for the component Uk, the 
complementary subset §(k) = {1, . . . M}\E(k) is called the set of static direc- 
tions for the component u^. The dependent variables Wfc(£) = Wfc(£i> • • • , Cm) 
are thought of as attached to the cubes of dimension adjacent to the 

point £. 

The functions /| ■ are defined on open domains Dfcj(e) C X and are of 
class C°° there. We suppose that the following conditions are satisfied. 

(F) The domains Ofcj(e) for 1 < j < m (that is for all directions which 
become continuous in the limit e — > 0) exhaust X as e — > 0, i.e. any 
compact subset K C X is contained in Dfcj(e) for all k, for all j G £(&), 
j < m, and for all e small enough. Similarly, there exists an open set 
Del such that the domains Dfcj(e) for m + 1 < j < M corresponding 
to the directions remaining discrete in the limit, exhaust D as e — > 0. 
The functions locally converge to the respective in C e (X), resp. 
in C e (B>), for any £, with the rate of convergence O(e). 

The existence of the subset D reflects the fact that transformations of conju- 
gate, resp. orthogonal, systems are not defined for all initial data. The set D 
consists of "good" initial values, where the transformations are well defined. 

Goursat Problem 1 Given functions Uf, on 7"§^) subspace spanned by 
the static directions of Uk), and r > 0, find a function u e : B e (r) — ► X 
that solves the system hll\) and satisfies the initial conditions u\ = on 

n (fc) ni? £ (r). 

In order for this problem to admit a unique solution for e > 0, this 
has to be the case just for one elementary M-dimensional cube B e (e). The 
corresponding condition is called consistency, and it turns out to be necessary 
and sufficient also for the global (formal) solvability of the above Goursat 
problem. The consistency condition can be expressed as 

Sj(SiUk) = 5i(8jUk), (12) 
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for any choice ofi^j from the respective £&, if one writes out all the involved 
5ii's according to the equations (JTTJ). The first step in writing this out reads: 

Sj(fkA u )) = S i(fk,j(u)), (13) 
In order for the left-hand side (say) of this equation to be well-defined, the 
function f^i is allowed to depend only on those ui for which j G &£. This has 
to hold for all j G £&, j ^ i, so we come to the condition that f^i depends 
only on those ue for which £& \ {i} C £^. If this is satisfied,then the final 
form of the consistency condition reads: 

ej 1 (jkA u + e jfj( u )) - fkA u )) = ^ [fkj(u + eifi(u)) - fk,j{ u )) > ( 14 ) 

where we use the abbreviation fi(u) for the .fT-vector whose kth component 
is fk,i( u )i if ^ G £&, and is not defined otherwise. 

It is readily seen that, if one prescribes arbitrarily the values of all vari- 
ables Uk at £ = (recall, Wfc(O) is actually associated with the respective 
static facet of dimension #(§&) adjacent to the point £ = 0), then the system 
defines uniquely the values for all variables on the remaining facets of 
the unit cube if and only if it is consistent. The following proposition shows 
that the local consistency assures also the global solvability of the Goursat 
problem; the only obstruction could appear if the solution would hit the set 
where the functions fkj are not defined. 

Proposition 1 Let e > ; and assume that the system Ml\) is consistent. 
Then either the Goursat problemU\has a unique solution u e on the whole of 
B e {r), or there is a subdomain Q C B e (r) such that the Goursat problem^ 
has a unique solution u e on Q, and u e (l;) ^ ©(e) at some point ( G fi. 

Finally, we turn to the whole family (0 < e < eo) of Goursat problems for 
equations with the respective data U^. The following result has been 
proven in our recent paper [B MSj . 

Theorem 2 Let a family of Goursat problems^ be given. Suppose that they 
are consistent for all e > 0, and satisfy the condition (F). 

1. Let the functions Ui,...,U^ be smooth, and, if m' > 0, assume that 
the point (6^(0), . . . , U^-(0)) belongs to the subset D. Then there exists 
r > such that the Goursat problem^ at e = has a unique solution 
u : B(r) X. 

2. If the Goursat data locally 0(e) -converge in C°° to the functions 
U®, then there exists e\ G (0,eo) such that the Goursat problems^ are 
solvable on B £ (r) for all < e < e\, and the solutions u e : B e (r) — > X 
converge to u : B(r) — * X with the rate 0(e) in C°°(B(r)). 
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3 Conjugate nets 



We give definitions of continuous and discrete conjugate nets and their trans- 
formations, formulate corresponding Goursat problems, and prove conver- 
gence of discrete conjugate nets towards continuous ones. 

3.1 Basic definitions 

Recall that M = m + m'. In Definitional below it is supposed that m' = 0, 
M = m, while in Definition El it is supposed that m' = l, M = m + 1. 

Definition 2 A map x : B (r) i-> R N is called an m- dimensional conjugate 
net in ~R N , if didjX G span(<9jX, djx) at any point £ G B Q (r) for all pairs 
1 < i ^ j < m, i.e. if there exist functions Cij, Cji : B (r) t— > K such that 



Definition 3 A pair of m- dimensional conjugate nets x,x + : Bo(r) > Mr is 
called a Jonas pair, if three vectors diX, diX + and Smx = x + — x are coplanar 
at any point £ G B (r) for all 1 < % < m, i.e. if there exist functions 
CiM, cm% '■ B (r) — > R such that 



Remarks. 

• x + is also called a Jonas transformation of x. One can iterate these 
transformations and obtain sequences (x^°\x^\ . . . , x^) of conjugate 
nets x^ : B (r) — ► M. N . It is natural to think of such sequences as of 
conjugate systems with m continuous and one discrete direction. We 
will see immediately that the notion of fully discrete conjugate systems 
puts all directions on an equal footing. 

• A Combescure transformation x + of a; is a Jonas transformation for 
which vectors dix + and diX are parallel, i = l,...,m. This class of 
transformations is singled out by requiring c^,/ = in equation (fTB]). 

Discrete conjugate nets were introduced by Doliwa and Santini DS1 . 

Definition 4 A map x : B e (r) i— > ~R N is called an M -dimensional discrete 
conjugate net in M. N , if the four points x, TiX, TjX, and TiTjX are coplanar 
at any £ G B e {r) for all pairs 1 < i ^ j < M , i.e. if there exist functions 
c^ : B e {r) — ► R such that 



X . 



(15) 



diX + - diX = CMidiX + c iM (x + - x) . 



(16) 



(17) 
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If m! — 1, M — m + 1, then the M -dimensional discrete conjugate net, 
considered as a pair of functions x,x + : Bq{t) \— > M N , is also called a Jonas 
pair of m- dimensional discrete conjugate nets. 

3.2 Hyperbolic equations for conjugate nets 

Introducing M new functions Wi : B e (r) — > M N , we can rewrite (j!7j) as system 
of first order: 

6iX = Wi, (18) 
5iWj = cjiWi + CijWj , i=£ j, (19) 
SiC k j = {TjC ik )c kj + {TjC kl )cij - {nCk^Cij, i^ j ^k^i. (20) 

For a given discrete conjugate net x : B e (r) Mr, eq. (|18jl defines the 
functions Wj, then eq. (fTHj) reflects the property (|T7J) . and eq. pOjl is just a 
transcription of the compatibility condition Si(SjWk) = 5j{5iW k )- Conversely, 
for any solution of (JT%j) -()20| ) the map x : B e (r) — > IR has the defining 
property (fTTj) and thus is a discrete conjugate net. So, discrete conjugate 
nets are in a one-to-one correspondence to solutions of the system (|18 |) -(|20 j) . 

This system almost suits the framework of Sect. 12. 2| the only obstruction 
is the implicit nature of the equations (|2T)|) (their right-hand sides depend 
on the shifted variables like TiC k j which is not allowed in (fTTj) . We return 
to this point later, and for a moment we handle the system ()18|) -(I20 J) just 
as if it would belong to the class (jlljl . In this context, we have to assign to 
every variable a Banach space and static/evolution directions. Abbreviating 
w = (w±, . . . , wm) and c = (c^)^, we set 

X = R N {x} x (M N ) M {w} x M M ( M ~V{c}, 

with the norm from ()10|). We assign to x no static directions, to the only 
static direction i, and to c^ two static directions i and j. The corresponding 
Goursat problem is now formulated as follows. 

Goursat Problem 2 (for discrete conjugate nets). Given: 

• a point X £ E M N , 

• M functions W[ : D 5 - i— > M N on the respective coordinate axes 7], and 

• M(M — 1) functions C\a ■ O 5 ^- i— > M on the coordinate planes 3^-, 

find a solution (x e , w e , c e ) : B e (r) i— > X to i/ie equations < f7^)) -< T^) satisfying 
the initial conditions 

x*(0)=X*, wt\ n =Wt, r,iv : C], (21) 
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We discuss now the consistency of the discrete hyperbolic system (j!8|) - 
([20)1 . First of all, note that for any triple of pairwise different indices k) 
the equations of this system involving these indices only form a closed sub- 
set. In this sense, the system (IT£ j) - (|2T)|) consists of three-dimensional building 
blocks. One says also that it is essentially three-dimensional. Correspond- 
ingly, in the case M = 3 consistency is not an issue: all conditions to be 
verified are automatically taken into account by the very construction of the 
equations ifTSjl - ljSnjl . 

In this case M = 3 the initial data for an elementary cube point 
x(0), three vectors u>i(0), u?2(0), ^(0) which are thought of as attached 
to the edges incident to the point £ = and parallel to the correspondent 
axes, and six numbers %■(()), 1 < i ^ j < 3, which are thought of as 
attached to the plaquettes incident to the point £ = and parallel to the 
coordinate planes T^. From these data, one constructs first with the help 
of (fTHjl the three points (0) + ejWj(O), then one calculates TiWj for 

i 7^ j with the help of (fHJj). and further TjT^x by means of (fTHj). Notice 
that the consistency assured by the symmetry of the right- 

hand side of (|T^|) with respect to the indices i and j. Thus, the first two 
equations (fT%|) . (fEU) allow us to determine the values of x in seven vertices 
of an elementary cube. The remaining 8-th one, t^t^x, is now determined 
as follows. Calculate TiCkj from the (linearly implicit) equations (|2T)|) . then 
calculate TjTiWk from (fTUj) . and then tit 2 t 3 x from (JTHJ). Again, the consistency 
requirement Tj^Wk = TiTjWk is guaranteed, per construction, by eqs. (ffilj) . 
Geometrically, the need to solve the linearly implicit equations for T^Ckj is 
interpreted as follows. The point t^t^x is determined by the conditions 
that it lies in three planes TiUjk, 1 < % < 3, where Tjlljfc is defined as the 
plane passing through three points (rjX, TjTjX, r fc Tjx), or, equivalently, passing 
through TjX and spanned by TjiUj, TiWk- All three planes r 1 n 2 3, T 2 n 13 and 
r 3 rii2 belong to the three-dimensional space through x(0) spanned by Wi(0), 
1 < i < 3, and therefore generically they intersect at exactly one point. The 
above mentioned linearly implicit system encodes just finding the intersection 
point of three planes. 

Turning now to the case M > 3, we face a non-trivial consistency prob- 
lem. Indeed, the point r^r^Tex can be constructed as the intersection point 
of three planes TiTjllke, TiTkHje, TiTiUjk, where the plane TiTjTlki is defined 
as the plane passing through three points {t%t jX, TiTjT^x, t$t jT^x) , or, equiv- 
alently, passing through T^TjX and spanned by TiTjWk, TiTjW£. So, there are 
four alternative ways to determine the point Tj Tj t^tix, depending on which 
index plays the role of i. 

Proposition 2 The system (Q3)-(d?J) is consistent for M > 3. 
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Proof. The geometric meaning of this statement is that the above mentioned 
four ways to obtain t^tjt ^Tpx lead to one and the same result. This can be 
proved by symbolic manipulations with the equations (fT8|) - (j20|) . but a geo- 
metric proof first given in |Dglj is much more transparent. In the construc- 
tion above, it is easy to understand that the plane TiTjUke is the intersection 
of two three-dimensional subspaces TjlL,^ and TjU^e in the four-dimensional 
space through x(0) spanned by u>i(0), Wj(0), Wk(0), u^(0), where the sub- 
space TiUjki is the one through the four points (r ix, t tTjX, t it ^x, t%t ix) , or, 
equivalently, the one through TjX spanned by TiWj, TiWk and TiW£. Now the 
point TiTjTkTfX can be alternatively described as the unique intersection point 
of the four three-dimensional subspaces TjlL,^, TjUi^e, T^U^i and r^ITjjfc of 
one and the same four-dimensional space. ■ 

Now, we turn to the above mentioned feature of the system (fT%|) - (f2T)|) . 
namely to its implicit nature. The geometric background of this feature is 
the following. Though there is, in general, a unique way to construct the 
8-th vertex t^t^x of an elementary hexahedron out of the the seven vertices 
x, TiX and TiTjX, or, equivalently, out of x, Wi, Cij, it can happen that the 
resulting hexahedron is degenerate. This happens, e.g., if the point r^jx 
lies in the plane TkHij through (TkX,TiTkX,TjTkx). This is illustrated on Fig. 

While the example on the left of Fig. 0] is fine, the right one does not 




Figure 4: A non-degenerate and a degenerate hexahedron. 

have the combinatorics of a 3-cube (the top-right edge has degenerated to 
a point). As a result, some of the quantities TjC^- in equation ()20|) are not 
well-defined; this is reflected by the fact that the equations pUj) are implicit 
and have to be solved for TiCkj, which is possible not for all values of c but 
rather for those belonging to a certain open subset. 
To investigate this point, rewrite ()20|) as 

SiC kj = F( ijk) (c) + < jidjC, i.-)c I,;, + < Mj*'ki)c,j - < i{<\<'kj)c,j (22) 
Introducing vectors 5c and F(c) with MiM — 1)(M — 2) components labelled 
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by triples (ijk) of pairwise distinct numbers 1 < i,j,k < M, 





we present the above equations as 



(l-Q^c))6c = F(c), 



(23) 



with a suitable matrix Q e (c). We need to find conditions for 1 — Q e (c) to be 
invertible. 

3.2.1 One conjugate net 

If ml = 0, then all entries of Q e (c) are either zero or of the form ecjj. 
Therefore, on any compact K C X we have: ||0/(c)|| < Ct with some 
C = C(K) > 0, hence 1 — Q e (c) is invertible for e < eo = I/O, and the the 
inverse matrix (1 — 0/(c)) _1 is a smooth function on K and 0(e)-convergent 
in 0°°(K) to 1. So, the system (J2"3~j) is solvable, and 



For any solution (x, w, c) : Bq(t) — > X of this system the function x : Bq(t) — > 
is a conjugate net, because equations (|2*lj) and (|23|) yield the defining 
property (fTHj) . Conversely, if a; : -Bo( r ) — > is a conjugate net, then equa- 
tion ()24j) defines the vectors u>j, these fulfill (|2*3jl since this is the defining prop- 
erty of a conjugate net, and (j2T?j) results from equating di(djWk) = dj(diWk). 
So, the following is demonstrated: 

Lemma 1 If w! — 0, then the Goursat problem [H for discrete conjugate 
nets fulfills the condition (F). The limiting system for e = 0, consisting of 
\24\) -ft£ty)- describes continuous conjugate nets. 



6~iC kj = F iijk) (c) + 0(e) = c ik c kj + c ki Cij - c kj c i:j + 0(e) 



where the constant in 0(e) is uniform on any compact set. 
The limiting equations for (|18 )) -([20 )) in this case are: 



diX 
diWj 

9i C k j 




(24) 
(25) 
(26) 
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3.2.2 Jonas pair of conjugate nets 

In the case m' = 1 the matrix 1— Q e (c) is no longer a small perturbation of the 
identity, since eu = 1- More precisely, the matrix Q e (c) is block-diagonal, its 
(^) diagonal blocks Q\ij k } (c) of size 6x6 correspond to non-ordered triples 
of pairwise distinct indices {i,j, k}, with rows and columns labeled by the 
six possible permutations of these indices. Blocks for which k ^ M are of 
the form considered before - their entries are either zero or of the form ec^. 
Blocks where, say, k = M, admit a decomposition Q e ujM\{c) = A(c) + eB(c), 
where the matrices A and B do not depend on e. It is not difficult to calculate 
that det(l — A(c)) = (1 + CMi)(l + CMj)- So 1 — <5{ ijM }( c ) * s invertible if 
CMi,CMj 7^ —1? provided e is small enough. The inverse is smooth in a 
neighborhood any such point c, and 0(e)-convergent to (1 — A(c)) -1 . 
As formal limit of (fTSj) and (|T§]) . we obtain (i ^ j): 

diX = Wi (1 < i < m), (27) 

S M x = w M , (28) 

diWj = CjiWi + CijWj (1 < i < m), (29) 

5 M Wj = c jM w M + c M jWj. (30) 

There are also four different limits of equation ()2U|). depending on which 
directions i and j are kept discrete. We shall not write them down; they 
are easily reconstructed as the compatibility conditions for (}2"9*j) . (jHOJ), like 
SM(diWk) = di(6MWk) etc. Comparing (f2"Tj) - (j30|) with the definition of the 
Jonas transformation (fTT)|) . we see that the following is demonstrated. 

Lemma 2 Ifm' = 1, then the Goursat problem^ for discrete conjugate nets 
also fulfills the condition (F), with the domain 

D = {(x,w,c) G X | c Mi 7^ —1 for 1 < i < m}. 

The limiting system at e = describes Jonas pairs of continuous conjugate 
nets. 

3.3 Approximation theorems for conjugate nets 

We are now ready to formulate and prove the main results of this chapter. 

Theorem 3 (approximation of a conjugate net). Let there be given: 

• m smooth curves Xi : iPj i— > R , i = 1, . . . ,m, intersecting at a common 
point X = Xi(0) = • • • = X m (0); 
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• for each pair 1 < i < j < m, two smooth functions Cij, Cji : D 5 ^- t— > R. 
Then, for some r > 0: 

1. There is a unique conjugate net x : B (r) — > R such that 

x\ % =Xi, Cijfo^Cij. (31) 

2. The family of discrete conjugate nets {x e : B^(r) — > M. N } 0<e<£1 uniquely 
determined by requiring 

x*h t =X t \ T ., c%\ n = C i5 \ n , , (32) 

where c\j are the rotation coefficients of the net x € , 0(e) -converges in 
C°°(B (r)) to x. 

Proof. It is easy to reformulate the data (|31j) . (f3*2"|) into valid initial data 
for the Goursat problem on X for the systems (jHJ)-(j2Sl), resp. ([TH jl -pU jl . 
Namely, at e = one translates ()31|) into the equivalent data 

x(0)=X, Wi\y = diXi, r,/ :) ., (',,. (33) 

while for e > one translates ()32|) into the equivalent data 

x* (0) = X , < \ n = 8iXi \ n , 4 \ n , = C l3 \ n , . (34) 

Now the theorem follows directly from Lemma Q and Theorem 121 ■ 

Theorem 4 (approximation of a Jonas pair). Let, in addition to the 
data listed in Theorem^ there be given: 

• m smooth curves X^ : CP, — > M. N , i = l,...,m, such that all Xf 
intersect at a common point X + = X± (0) = . . . = X+(0) ; and such 
that for any 1 < i < m and for any point £ e CPj the three vectors diXi, 
diX^ and S^fXi = Xf — Xi are coplanar. 

Assume that X + — X is not parallel to any of the m vectors diXf(0) . Define 
the functions Cmu CiM '■ CP« t— > R by the formula 

diX+ - d t X t = C Mi diXi + C m {Xt - X t ). (35) 

Then, for some r > 0: 

1. In addition to the conjugate net x : B (r) — ► R^ defined in Theorem^ 
there is its unique Jonas transformation x + : Bq[t) — > M. N such that 

x + W = X+ (l<i<m); (36) 
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2. The family {x e : B e (r) —>■ M. N } 0<t<tl of discrete M -dimensional con- 
jugate nets (M = m + 1), uniquely determined by the requirements 

x e \y t = Xi \ n , 4 \ n j = Cij \ n . , (l<i^j<m) (37) 

and 

x e (0, 1) = X + , c e Mi \ n = C M i\vf, c\ M \ n = C lM \?i (l<i<m) 

(38) 

0(e)- converges in C°°(B(r)) to the Jonas pair x,x + , in the sense that 
x e (-, 0) — > x, x e (-, 1) — > x + 

Proof. Consider the Goursat problem for M-dimensional discrete conjugate 
nets with the initial data 

x € (0)=X, wt\ n = SiXify, w e M (0) = X + -X, (i = l,...,m) (39) 

and 

(40) 

The statement of the theorem follows by applying Theorem El to this situa- 
tion. This, in turn, is possible due to Lemma|2]and the following observation. 
The condition of the theorem yields that (£) — JQ(£) is not parallel to any 
of diXi~(£) not only at £ = 0, but also in some neighborhoods of zero on 
the corresponding axes JV Now one deduces from the definition (|35j) of the 
rotation coefficients CW, Cm% that in the same neighborhoods Cm% ^ — 1- 
Hence, the data of our Goursat problem belong to the set D of Lemma EJ ■ 

Letting not one but two or three directions of the conjugate net remain 
discrete in the continuous limit (so that ml = 2 or ml = 3), one arrives at 
the permutability properties of the Jonas transformations. 

Theorem 5 (permutability of Jonas transformations). 

1. Given an m-dimensional conjugate net x(-,0,0) : B Q (r) \—> M. n and 
its two Jonas transformations x(-,l,0) : B Q (r) i— > M, N and x(-,0,l) : 
A)( r ) l— ► K^, there exists a two-parameter family of conjugate nets 
x(-, 1, 1) : Bo(r) h-» ~R n that are Jonas transformations of both x(-, 1, 0) 
andx(-, 0, 1). Corresponding points of the four conjugate nets are copla- 
nar. 
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2. Given three Jonas transformations 

X{; 1, 0, 0), x{; 0, 1, 0), x{; 0, 0, 1) : B {r) h+ R N 

of a given m -dimensional conjugate net x(-, 0,0,0) : B (r) i— > R^, as 
well as three further conjugate nets 



x(-,l,l,0),x(-,0,l,l),x(-,l,0,l): B (r) 



stzc/i i/iai 1, 1, 0) is a Jonas transformation of both x(-, 1, 0, 0) and 
x(-, 0,1,0) etc., there exists generically a unique conjugate net 



x{-,l,l,l) :B (r) 



which is a Jonas transformation of all three x(-, 1, 1, 0), x(-, 0, 1, 1) and 
£(•,1,0,1). 

These results are proved in the exactly same manner as Theorem EJ if one 
takes into account the mult i- dimensional consistency of discrete conjugate 
nets (Proposition |2J). 



4 Orthogonal systems 

We give the definitions of orthogonal systems and their transformations, for- 
mulate corresponding Goursat problems, and prove convergence of discrete 
orthogonal systems to continuous ones. 

4.1 Basic definitions 

Definition 5 A conjugate net x : Bq (r) ^ R N is called an m-dimensional 
orthogonal system in ~R N ( an orthogonal coordinate system if m = N), if 

diX ■ djX = 0. (41) 

Definition 6 A pair of m-dimensional orthogonal systems x,x + : Bo(r) i— > 
W N is called a Ribaucour pair, if for any 1 < i < in the corresponding 
coordinate lines of x and x + envelope a one-parameter family of circles, i.e. 
if it is a Jonas pair of conjugate nets and 

<l '~ + ^f\.( x +- x )=Q. (42) 



diX + \ \dix\ 



Remarks. 
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• x + is also called a Ribaucour transformation of x; Ribaucour trans- 
formations can be iterated, and our results generalize to any finite 
sequence of such transformations (cf. the remarks after Definition EI). 

• Two-dimensional (m = 2) orthogonal systems are called C-surfaces. 
C-surfaces in M. 3 are characterized as surfaces parametrized along cur- 
vature lines. 

• There are different definitions of Ribaucour transformations in the lit- 
erature (cf. [H ). We have chosen the one that is best suited for our 
present purposes. 

Definition 7 A map x:B e (r) t— > M. N is called an M- dimensional discrete or- 
thogonal system inM. N , if the four points x , TiX, TjX andr^x are concircular 
for all £ E B e (r) and all 1 < j ^ j < M. 

If m' = 1, M = m + 1, then the M -dimensional discrete orthogonal 
system, considered as a pair of functions x,x + : Bq{t) \— > M. N , is called a 
Ribaucour pair of m- dimensional discrete orthogonal systems. 

Remark. 

• Also in the discrete case we shall use the term C-surfaces for two- 
dimensional (m = 2) orthogonal systems. 

So, orthogonal systems form a subclass of conjugate nets subject to a certain 
additional condition (orthogonality of coordinate lines, resp. circularity of 
elementary quadrilaterals). 

The classical description of continuous orthogonal systems is this (in the 
following equations it is assumed that i ^ j ^ k ^ i): 



diX = h^, (43) 

diVj = f3jiVi, (44) 

dihj = hiPij, (45) 

dipkj = PkiPij, (46) 

diPij + djPji = -diVi ■ djVj . (47) 



Eq. (|4*Hj) defines a system {vi}i<i< m of m orthonormal vectors at each point 
£ G B (r); the quantities hi = \dix\ are the respective metric coefficients. 
Conjugacy of the net x and normalization of t>j imply that eq. (}4*4"j) holds 
with some real-valued functions /3y. Eq. (|43|) expresses the consistency 
condition diidjx) = dj(dix) for i ^ j. Analogously, eq. (|4T)j) expresses the 
consistency condition di(djVk) = dj(diVk) for i ^ j ^ k ^ i. The m equations 
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and (™) equations (}4l)j) are called the Darboux system; they constitute 
a description of conjugate nets alternative to the system (}2"4"|) - (f2l)j) we used 
in Sect. El The orthogonality constraint is expressed by (™) additional 
equations (|4*7|) . derived from the identity didj(vi, Vj) = 0. In the case m = N 
the scalar product on the right-hand side of (jUj) can be expressed in terms 
of rotation coefficients /3y only: 

+ djPji + PkiPkj = 0. (48) 

The Darboux system (|45jl . (|4fi|) together with the eqs. (|4"%|) forms the Lame 
system. 

The classical approach as presented is based on the Euclidean geometry. 
However, the invariance group of orthogonal systems is the Mobius group, 
which acts on the compactification ~R N U {oo} ~ S N , rather than on M. N . 
This is a motivation to consider orthogonal systems in the sphere and to 
study their Mobius-invariant description. In particular, this enables one to 
give a frame description of orthogonal systems which can be generalized to 
the discrete context in a straightforward manner. This turns out to be the 
key to derivation of the discrete analogue of the Lame system. 



4.2 Mobius geometry 

We give a brief presentation of Mobius geometry, tailored for our current 
needs. A more profound introduction may be found in |Hj . 

The iV-dimensional Mobius geometry is associated with the unit sphere 
C 1R 7V+1 . Fix two antipodal points on E> N , po = (0, . . . , 0, 1) and p^ = 
—po- The standard stereographic projection a from the point is a con- 
formal bijection from §f = \ {p^} to M that maps spheres (of any 
dimension) in to spheres or affine subspaces in Mr. Its inverse map 



_ i — 1 I x I 

1 + \xr 1 + \x 12 



2 



lifts an orthogonal system in R , continuous or discrete, to an orthogonal 
system in S^. 

The group M(iV) of iV-dimensional Mobius transformations consists of 
those bijective maps fi : — > that map any sphere in to a sphere of 
the same dimension; it then follows that \x is also conformal. 

Fact 1 The notion of orthogonal system on §> N is invariant under Mobius 
transformations. 
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Mobius transformations can be embedded into a matrix group, namely the 
group of pseudo-orthogonal transformations on the Minkowski space R Ar+1 ' 1 , 
i.e. the (N + 2)-dimensional space spanned by {ei, . . . , ejv +2 } and equipped 
with the Lorentz scalar product 

N+2 N+2 N+l 

^ Ui^i, ^2 w i e i) = UkVk ~ U N+2VN+2- 
i=l j=l ' k=l 

To do this, a model is used where points of are identified with lines on a 
light cone 

L N+1 = {u e R N+1>1 I (u, u) = 0} C R N+1 '\ 
This identification is achieved via the projection map 

7T : R^ 1 ' 1 \ (R N+1 x {0}) i — > R N+1 , 

(u 1 ,...,U N+1 ,U N+2 ) ^ («l/«JV+2 , • • • ,U N+ i/u N+2 ), 

which sends lines on £L N+1 to points of S . In particular, the lines through 

e o = ^( e N+2 + ^n+i) and = ~(e N+2 - e N+1 ) 

are mapped to the points po and p^, respectively. 

The group + (N + 1,1) of genuine Lorentz transformations consists of 
pseudo-orthogonal linear maps L which preserve "the direction of time" : 

(L(u), L(v)) = (u, v) V«, v G R N+1 >\ (L(e N+2 ), e N+2 ) < 0. (49) 

For L G + (N + 1,1), the projection n induces a Mobius transformation of 
via fi(L) = 7i o L o 7T _1 . Indeed, L is linear and preserves the light cone 
because of (J4*9~j) . so fi(L) maps to itself. Further, L maps linear planes 
to linear planes; under tt, linear planes in R^" 1 " 1 ' 1 project to afline planes in 
R^" 1 " 1 . Since any sphere in is uniquely represented as the intersection of 
an afline plane in M. N+1 with C R^"^ 1 , we conclude that spheres of any 
dimension are mapped by fi(L) to spheres of the same dimension, which is 
the defining property of Mobius transformations. It can be shown that the 
correspondence L \— > (jl(L) between L G + (N + 1, 1) and M,(N) is indeed 
one-to-one. 

Recall that we are interested in orthogonal systems in R^; their lifts 
to actually do not contain p^. Therefore, we focus on those Mobius 
transformations which preserve p^, thus corresponding to Euclidean motions 
and homotheties in R^ = cr(S^). For these, we can restrict our attention to 
the section 

% = |u G R^ 1 ' 1 : (it, u) = 0, (u, eoo ) = -1/2} = L N+1 n (e + e^). 
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(We use here the notation u H for the hyperplane orthogonal to u.) The 
canonical lift 

A : R N i — > X, 

x h-> x + e + l^peoo, 

factors the (inverse) stereographic projection, a -1 = 7roA, and is an isometry 
in the following sense: for any four points X\, . . . , X4 G M. N one has: 

(X(xi) - X(x 2 ), \{x 3 ) - A(x 4 )) = (xi - x 2 ) ■ (x 3 - x A ). (50) 

We summarize these results in a diagram 

% ^ c m jv+i,i 

At i* 

Denote by 0+ (iV + 1,1) the subgroup of Lorentz transformations fixing the 
vector (and consequently, preserving %). The same arguments as before 
allow one to identify this subgroup with the group E(N) of Euclidean motions 
ofR* 

Fact 2 Mobius transformations of S> N are in a one-to-one correspondence 
with the genuine Lorentz transformations ofR N+1,1 . Euclidean motions of 
W N are in a one-to-one correspondence with the Lorentz transformations that 
fix . 

Our next technical device will be Clifford algebras which are very convenient 
to describe Lorentz transformations in R^ 4 " 1 ' 1 , and hence Mobius transfor- 
mations in 8 . Recall that the Clifford algebra C(iV + l, 1) is an algebra over 
R with generators ei, . . . , eN+2 G R^" 1 " 1 ' 1 subject to the relations 

uv + vu = -2(u,v)l = -2(u,v) V«,dgE w . (51) 

Relation (|51jl implies that u 2 = —(u,u), so any vector u G R^ 4-1,1 \ L N+l 
has an inverse w _1 = —u/(u,u). The multiplicative group generated by 
the invertible vectors is called the Clifford group. It contains the subgroup 
Pin(N + 1,1) which is a universal cover of the group of Mobius transforma- 
tions. We shall need the genuine Pin group: 

<K = Pin + (N + 1, 1) = {«!•• • u n I u\ = -1}, 

and its subgroup generated by vectors orthogonal to eoo: 

^00 = {ux---u n \ u\ = -1, (ui, e^) = 0} C %. 
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"K and Oi^ are Lie groups with Lie algebras 



fj = spin(N + 1,1) = span^e^ : i,j G {0, 1, . . . , N, oo}, i ^ jj, (52) 
f)oo = spin^N + 1, 1) = spaniel : i, j G {1, . . . , N, oo}, % ^ j }. (53) 



In fact, % and "K^ are universal covers of the previously defined Lorentz 
subgroups + (N + 1,1) and 0+(iV + 1,1), respectively. The covers are 
double, since the the lift ip{L) of a Lorentz transformation L is defined up to 
a sign. To show this, consider the coadjoint action of ip G !K on v G ]R 7V+1,1 : 



Thus, A u is, up to sign, the reflection in the (Minkowski) hyperplane w H 
orthogonal to u; these reflections generate the genuine Lorentz group. The 
induced Mobius transformation on S w (for which the minus sign is irrelevant) 
is the inversion of 8^ in the hypersphere 7r(-u H n £j N+1 ) C S n ; inversions in 
hyperspheres generate the Mobius group. In particular, if u is orthogonal 
to e^, then A u fixes e^, hence leaves % invariant and induces a Euclidean 
motion on ~R N , namely the reflection in an afline hyperplane; such reflections 
generate the Euclidean group. 

Fact 3 Mobius transformations ofE> N are in a one-to-one correspondence to 
elements of the group CK/{±1}; Euclidean transformations of ~R N are in a 
one-to-one correspondence to elements of < K tXl /{±.\}. 

4.3 Orthogonal systems in Mobius geometry 

We use elements ip G 'Koo as frames when describing orthogonal systems in 
the Mobius picture. The vectors A^(ej) for % — 1, . . . , N form an orthogonal 
basis in T% at the point A^(e ) G %. For an orthogonal system x e : B e (r) i— > 
%, continuous (e = 0) or discrete (e > 0), with [m! > 0) or without (m! = 0) 
transformations, consider its lift to the conic section %: 



(54) 



Obviously, for a vector u with u 2 



— 1 one has: 




(55) 



x' = \ox' : B e (r) — > X. 



(56) 
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4.3.1 Continuous orthogonal systems 

Define Vi : Bq(t) — > T3C for z = 1, . . . , m by 

dd = hiVi, hi = \dix\. (57) 

The vectors v% are pairwise orthogonal, orthogonal to e^, and can be written 

as 

Vi = Vi + 2(x,v i )e 00 (58) 

with the vector fields Vi from ([43)1 . As an immediate consequence, these 
vectors satisfy 

with the same rotation coefficients (3ji as in (|44j) . 

Definition 8 a) Given a point x G % and vectors t)j G T^IJC, 1 < i < m, 
we call an element ip G Jioo suited to (x; Vi, . . . , v m ) if 

A/,(e ) = x, (60) 
Ai/,(e k ) = v k (1 < k < m), (61) 

b) A frame ip : B (r) t— > Ji^ is called adapted to the orthogonal system 
x, if it is suited to (x; V\, . . . , v m ) at any point £ G B (r), and 

diA^e k ) = f3 kt Vi (1 < k < N, 1 < i < m). (62) 

with some scalar functions f3 k i- 

Proposition 3 For a given continuous orthogonal system x : B (r) \— > %, 
an adapted frame ip always exists. It is unique, if m = N. If m < N , then 
an adapted frame is uniquely determined by its value at one point, say £ = 0. 
An adapted frame ip satisfies the following differential equations. 

• Static frame equations: 

diip = —e^Si (1 < i < m) (63) 

with Si = {l/2)diVi. 

• Moving frame equations 

= -Siei?p (1 < i < m) (64) 
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with 

where the functions and hi constitute a solution to the Lame system 
(in the following equations the indices 1 < i,j < m and 1 < k < N are 
pairwise distinct): 

dihj = hifaj, (66) 

diPhj = PkAj , (67) 

diPy + djpji - - ^ ■ (68) 

Conversely, if and hi are solutions to the equations \6b]) - KUB) . £/ien i/ie 
moving frame equations \64\) are compatible, and any solution is a frame 
adapted to an orthogonal system. 

Proof. For an orthogonal coordinate system (m = N) the adapted frame is 
uniquely determined at any point by the requirements fl60l) and (jolj) . while 
eqs. ((02)) follow from (jfilj) . If m < N, extend the m vectors Vi{0), 1 < i < m, 
by N — m vectors Vk(0), m < k < N, to a orthonormal basis of 
There exist unique vector fields v/. : B Q (r) i— > m < k < N with these 

prescribed values at £ = 0, normalized (-0^ = —1), orthogonal to each other, 
to 0j, 1 < i < m, to x and to e^, and satisfying differential equations 

didk = PkiVi , 1 < i < m, m < k < N (69) 

with some scalar functions (3^- Indeed, orthogonality conditions yield that 
necessarily (3 ki = — (SjWj, Eqs. (JBTJj) with these expressions for f3 ki form a 
well-defined linear system of first order partial differential equations for dk- 
The compatibility condition for this system, dj(diVk) = dj(diVk) (1 < i ^ 
j < m), is easily verified using (JHHJ). The unique frame if) : B (r) i— > "H^ with 
A/>( e o) = &i A/>( e fc) = Vk (1 < k < N) is adapted to x. 

Next, we prove the statements about the moving frame equation. The 
equations 

<9iA/,(e ) = hiA^(ei), diA,p(e k ) = ^A^e,) 
are equivalent to 

[e , {dii))i)~ 1 } = K&i , [e fc , (d^)^ -1 ] = Phi^i ■ 

Since, by (JH3J) . (diip)" 1 ^ is spanned by bivectors, we come to the represen- 
tation of this element in the form (diip)' 1 ^ = — S^e, with S% as in (|65jl . The 
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compatibility conditions for the moving frame equations, di(djip) = dj(diijj), 
are equivalent to 

OjSiGj -\- diSj&i + Si(&iS j&i) + Sj{&jSi&j) = 0, 

which, in turn, are equivalent to the system (ffiT)J) - (pTHJ) . Conversely, for a 
solution ij) : B (r) i— > "K^ of the moving frame equations, define x = A^(e ) 
and Vk = A^{e k ). Then diX = hidi, and orthogonality of t>j yields that x is 
an orthogonal system. 

We conclude by showing the static frame equation. If ip is adapted, then 
(|f)5jl in combination with the orthogonality relation of the i>, yields: 

Si = A e .^(Si) = - f3 ki A e -^(e k ) — hjAe^epq) 

k^i 

= --^2(3 k iA^(e k ) + hiA^,(e ao ) = —^PkiVk + ^ e oo = 7^^. 

The last equality follows from (3 ki = — (diVi,v k ) and the fact that the e^ 
component of diVi is equal to 2/i, (the latter follows from (JHHjl ). ■ 



4.3.2 Discrete orthogonal systems 

Now consider a discrete orthogonal system x : B e (r) 1— > 3C. Recall that the 
cases m/ = and m' > are treated simultaneously, so that 6j = e for 
1 < z < m, and 6j = 1 for m + l<i<m + m' = M. Appropriate discrete 
analogues of the metric coefficients hi and vectors £>j are 

hi = \5ix\, di = h^SiX. (70) 

These are unit vectors, i.e. vj = —1, representing the reflections taking x to 
TiX in the Mobius picture, cf. Fig. |SJ 

T { x = -Aq.x = x + eihiVi . (71) 

It is important to note that vectors v,i are not mutually orthogonal. Next, 
we derive equations that are discrete analogues of (fHHJ) . Since four points x, 
T{X, Tjx and TiTjx are coplanar, the vectors TjSiX and TiSjX lie in the span of 
5iX and SjX. The lift A is affine, therefore the same holds for t^'s: 

TiVj = {njiY 1 ■ (vj + pjiVi). (72) 

Here pji are discrete rotation coefficients, with expected limit behaviour p^ pa 

€i€jj3ji as e — ► 0, and 

n % = 1 + ^ + 2 Pji(vi, Vj) (73) 
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Figure 5: An elementary quadrilateral of a discrete orthogonal system 



is the normalizing factor that is expected to converge to 1 as e — > 0. Cir- 
cularity implies that the angle between 5iX and 5jX and the angle between 
TjSiX and TiSjX sum up to n: 

(TjVi, Vj) + (Vi, TiVj) = 0. (74) 

Under the coplanarity condition (f?2*j) the latter condition is equivalent to 

VjirjVi) + Viinvj) = 0. (75) 

Inserting (J72J) and (J73J) into (|74*jl . one finds that either (vi,Vj) = 1, i.e., the 
circle degenerates to a line (we exclude this from consideration), or 

Pij + Pji + 2(v u vj) = 0. (76) 

^From (|75|) and (|75|) there follows: 

n% = nji = 1 - PijPji. (77) 

So, the orthogonality constraint in the discrete case is expressed as the con- 
dition (|76| on the rotation coefficients; recall that in the continuous case this 
was a condition (jjTj), resp. (|38|). on i/ie derivatives of the rotation coefficients. 

The following definition mimics the static frame equations of the contin- 
uous case. 

Definition 9 A frame ip '■ B e {r) \— > [Kqo is called adapted to a discrete or- 
thogonal system x, if 

^(eo) = x, (78) 
= -e^ Vi (l<i<M). (79) 
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Existence of such a frame is a consequence of the circularity condition. In- 
deed, the consistency TjTiip = TiTjip is written as 

eiejipVj^TjVi) = eje^Vifavj). 

which is equivalent to ([75)1 . An adapted frame is uniquely defined by the 
choice of ipifS) (and thus is not unique, even if M = N). Indeed, eq. (IT9|) 
together with the definition (|7()j) of %i imply that if eq. (J7%|) holds at one 
point of B e (r), then it holds everywhere. 

Proposition 4 An adapted frame ip for a discrete orthogonal system x sat- 
isfies the moving frame equations: 



with 



where 



nip = (1 < i < M) (80) 



A e^(Vi) = N i e i + Tj^PkiGk ~ tiKeoo, 



6? 



and the functions flij and hi solve the following discrete Lame system (in the 
equations below 1 <i,j < M , 1 < k < N, and i ^ j ^ k ^ i): 

Tihj = [iij.i j) '•('.;/', • ( J i; I',,)- (83) 
TiPkj = (n-jitj)' 1 ■ (tjPkj + tiPkipij), (84) 
Tifaj = (i'j,< j) 1 • ('2\ii>,j ■ < j-hj). (85) 

Pij + Pji = NiPijej + NjPjiei - ^ ^ Pufikj- (86) 

where the abbreviation = (1 — PijPji) 1 ^ 2 is used, and are suitable real- 
valued functions. Conversely, given a solution and hi of the equations 
n£3\) - n?b}) . the system of moving frame equations (Q??J) is consistent, and its 
solution ip is an adapted frame of a discrete orthogonal system. 

Proof. From the definition E, = A~,(vi) it follows that E^ is a vector 
without e -component, i.e. admits a decomposition of the form (|%T|) . Again, 
from the definition of E and the moving frame equation (f?T)|) we find: 

A ei (Ej) = eiejipVjip^ejei , r^E,- = -e j e i V'^(Ti%)i)iV ; ~ 1 eie J - . 
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Upon using the identity Vi(vj + PjiVi)Vi = i)j + PijVi , which follows easily from 
(J7Bj), we can now represent eq. (|72|) in the following equivalent form: 

UjiiTiEj) = —A e . (Hj) - pijA^^i). (87) 

This equation is equivalent to eqs. (J53*j) - (J55)l . Next, the consistency Tjfai/j) = 
TiijjTp) of the moving frame equations (J80|) is equivalent to 

(TjE^CEi) + (rjEi)A ei (Ej) = 0. (88) 

Inserting the expressions for r^E,- and TjE from (|87j) into the left side give: 

A ei (S,)A ej (E) + A,(E)A,(£,) + pji(A ei (Ej)) 2 + Pij (Ae^i)) 2 

= NiP^i + NjPijej - (6^/2) ^ /3 ki (3 kj - (p„ + p^) = 0, 

the last equality being nothing but eq. (jSE}- Conversely, eqs. flHHj) - (J%o]) 
imply (|8Tj) . as well as the compatibility of the frame equations fjSOj) . So for an 
arbitrary initial value ^(0) G 3-Coo, the frame equations for ip : £? e (r) i— >■ < K co 
can be solved uniquely. Set x e = A^(eo), ■fij = A ei ^(Ej). ^From (|S7jl . (JHHJ) for 
the quantities E, there follow (|72|). (|74*j) for the quantities t>j. Using the fact 
that Ej has no e -component, one shows easily that there holds also (j71|) . 
Therefore, x is a discrete orthogonal system. ■ 

4.3.3 Ribaucour transformation of a continuous orthogonal 
system 

To describe a Ribaucour pair x, x + : Bq{t) \— > % of continuous orthogonal 
systems, one combines the descriptions of two previous subsections. Recall 
that in the present context M — m + 1. 

We denote by hf, vf, (3^ the corresponding objects for x + , defined as in 
(|H7j) . (JHnj) - Denote by % unit vectors such that 

x+ = -A iM (x), (89) 

cf. eq. (ITT]) . The defining property of Ribaucour transformations and 
the normalization of % imply: 

T> t + = (90) 
d^M = y(^ + + ^i) = oti (v{- (vM,Vi)vM), (91) 

with m auxiliary functions «j : -Bo( r ) ^ 
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Definition 10 A pair of frames : B (r) i— > zs called adapted to 

the Ribaucour pair x, x + , if ip is adapted to x and 

ip + = -e M ipv M . (92) 

Proposition 5 Let be a pair of frames adapted to the Ribaucour pair 

x, x + . Then the frame ip + is adapted to x + , and the following moving frame 
equations hold for ip + : 

= -Sfe^ (1 < i < to), (93) 
^ + = -S A fe M ^, (94) 



where 



= A~^(v M ) = N M e M + -z ^2 PkMe k - h M eoc, (96) 



2 

k^M 



N 2 M = 1 



\ E &m> (97) 



4 

k^M 



the functions hi, /3y and solve eqs. URfy - KUBj) . and the following system 

is satisfied (for 1 < i ^ j < m, 1 < k < N , and i ^ k ^ M ^ i): 

hf = hi + h M oti , (98) 
= pM+PkMVi, (99) 

Ptti = -PMi + 2N M a il (100) 
dih M = ^{h i + ht)i3 m , (101) 

diPkM = liPki + P^PiM, (102) 
diPiM = -^J2^ ki + (3 ki)PkM + N M (p Ml -Pt u ), (103) 

k^i,M 

diaj = \a i {p iS + ^). (104) 

Conversely, given solutions hi, and hf,/3^ of the system above with suit- 
able auxiliary functions on, the moving frame equations 4,9,^) . \94\) are com- 
patible, and the solution ip, ip + is a pair of frames adapted to a Ribaucour 
pair of orthogonal systems. 
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Proof. First we show that tp + is an adapted frame for x + . We have for 
1 < i < m, 1 < k < N: 

A/>+(e ) = ^ M A eM ^(e ) = -Ao M (x) = x + , 

A^+(ei) = A CM A eMl p(ei) = -A iM (vi) = vf , 

d i A il +{e k ) = di(Ai, M A eMl p(e k )) = ±d i (A i>M A i; (e k )) 

= ±di{A^{e k ) - 2{v M , A\, p (e k ))v M ) 

= ±((3ki - 2ai(v M ,A^(e k )))(vi - 2(v M , Vi)v M ) 

= ±(p k i-2ai(v M ,A^(e k )))v? = 0fcv+, 

where the minus sign applies iff k — M (and if k = 0, the latter calculation 
goes through almost literally, with replacing f3 ki , (3^ by hi, h + , respectively). 
Next, from Tim — ^eMtt(^) there follows: 

(3 kM = 2(E M ,e k ) = 2(v M ,A eM Tp(e k )) = -2(v M ,A^(e k )), 
h M = 2(S M ,e ) = 2(v M , A eM ^(e )) = -2(v M ,x), 
N M = (S M ,e A/ ) = (%, A eMl p(e M )) = (v M , A^(e M )). 

This proves eqs. (jUH|) - (jl00j) . Further, eqs. (jl(JH) - ()103J) are readily derived 
by calculating the i-th paartial derivative of the respective scalar product. 
Finally, eq. (jl()4j) comes from the consistency condition di(djVM) = ^(c^m)- 
Conversely, given a solution to the equations and (|HK|) - (jl()4j) . 

the moving frame equations are consistent, thus defining the frames tp, 
It follows from Proposition El that both x = A^(e ) and x + = A^+(e ) are 
orthogonal systems. Furthermore, eq. (JH2), yields the defining relations 
flHU-flHIJ) of the Ribaucour pair, with vm = A Bm ^(Em)- ■ 

5 Goursat problems and approximation for 
orthogonal systems 

It would be tempting to derive the theory of the Lame system ()43j) - (f4"Tj) 
from its discrete counterpart (|H3^) - (jHB|) . treating the latter as a hyperbolic 
system of the type (jllj) . However, it turns out that in dimensions M > 3 
this is hard to carry out, since one needs to enlarge the set of dependent 
variables and equations in a cumbersome manner. The way around is based 
on the following fundamental lemma, which allows one to take care of two- 
dimensional orthogonal systems (in the coordinate planes) only, and to use 
then the results for conjugate systems. 
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Lemma 3 a) If for a conjugate net x, continuous or discrete, its restriction 
to each plane DV,-, 1 < i ^ j < M, is a C-surface, then x is a (continuous or 
discrete) orthogonal system. 

b) If in a Jonas pair (x,x + ), continuous or discrete, the net x is an 
orthogonal system, and the corresponding coordinate curves x \j>. and x + \-j> t 
envelope one-dimensional families of circles, then (x,x + ) is a Ribaucour pair 
of orthogonal systems. 

Proof in the discrete case is based on the Miguel theorem, cf. Sect. HJ and 
can be found in |CDS| . For the continuous case it is a by-product of the 
proof of Theorem |U] below. ■ 

So, suppose that M — 2. Then the moving frame equations (j80|) and the 
system (|B3 jl - (|HRjl take the form (1 < i, j < 2, 1 < k < N, i ^ j ^ k ^ i): 

6iip = (— l - -)-^f3 ki e k e i + h i e 00 ei\'4) , (105) 

Sihj = ^hi + ^—^hi, (106) 

ejn €in 

8Aj = ^^-—Pu, (107) 

5i(3 kj = ^—^P kj + —Pki, (108) 
P12 + P21 = e 2 N l (3 12 + e 1 N 2 f3 21 - e^Q. (109) 
Here we use the abbreviations 

I' = '1 12 = "2j - a/1 - P12P21, 6 = 2 ^2PklPk2, Nf = 1 - j ^Pli . 

k>2 k^i 



(HO) 

We show how to re-formulate the above system in the hyperbolic form and 
to pose a Cauchy problem for it. 

5.1 Hyperbolic equations and approximation for 
C-surfaces 

First consider the case m = 2, m! = 0, ei = e 2 = e, corresponding to C- 
surfaces. Set 

P12 = eiV 1 /? 12 -( e 2 /2)(0- 7 ), (HI) 
P21 = eiV 2 /3 21 -( e 2 /2)(e + 7 ), (112) 
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with a suitable function 7 : y\ 2 1— ► R. It can be said that the auxiliary 
function 7 splits the constraint, making the Lame system hyperbolic. This 
splitting plays a crucial role for our approximation results. In the smooth 
limit 2 7 = <9i/?i 2 - d 2 /3 2 i (cf. eqs. ffTT4l. fTTK)l below). 
As a Banach space for the system take 

X = IK{^} x R^-^iP} x R 2 {h u h 2 } x M{ 7 } , 

where "K is the space of an (arbitrary) matrix representation of Spin(N+l, 1) 
(note that if ip G JCoo C 5i at some point, then eq. (|105j) guarantees that 
this is the case everywhere), (3 denotes the collection of fikj with k = 1,2, 
1 < j < iV", 7^ j. Obviously, 

iVi = l + 0(e 2 ), iV 2 = l + 0(e 2 ), n = l + 0(e 2 ), (113) 

where the constants in O-symbols are uniform on compact subsets of X. 

Both directions % — 1, 2 are assumed to be evolution directions for ■?/>; for 
Ttj and /3fcj there is one evolution direction j ^ i, while for the additional 
function 7 both directions i = 1,2 are static. 

Goursat Problem 3 (for discrete C-Surfaces). Given G "Koo, func- 
tions E\ : y\ R and Bj^ : 0* ^ R /or z = 1, 2, 1 < fc < iV ; jfe ^ i, and 
r e : T l2 h-> R, find a solution to the equations fffltyl - fllty) with 1777)) . 177^) . 

and with the Goursat data 

^(o) = * e , K\ n =Hi M n = B tn 7h 2 =r e . 

Lemma 4 T/ie hyperbolic system hlO^ - HUffy with Mll\) . hlltyl is consis- 
tent. Goursat problem^ for discrete C-surfaces satisfies condition (F). The 
limiting system for e = consists of equations ffifl) with (E2J), KUBj) . ^Ffi) for 

paiwise distinct indices 1 < i, j < 2, 1 < k < N, and 

dxfa = -(1/2) J2^k2 + 7, (H4) 

k>2 

d 2 (3 21 = -(l/2)^/3 fcl ^2-7, (H5) 

k>2 

which describe continuous C-surfaces. 

Proof. Consistency is easy to see: the only condition to be checked is 
d~id 2 ip = o~ 2 o~iip, but this has already been shown in Proposition |3J Also other 
statements are obvious from l| 113)1 . Notice that hyperbolic equations (j!14j) . 
(1115)) come to replace the non-hyperbolic orthogonality constraint (|fi8j) . ■ 
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Now we discuss the way to prescribe the data in the Goursat problem 
E] in order to get an approximation of a given smooth C-surface. Unlike in 
the case of general conjugate nets, it is not possible to prescribe the discrete 
curves X\ = x e \yt coinciding with their continuous counterparts X, = x \y i 
at the lattice points. However, it is still possible to achieve that X\ are 
completely determined by Xj. 

Given a curve Xj : CP^ i— > % ) one has the corresponding tangential vector 
field d{Xi, and therefore the function hi = |<9jXj| and the field of unit vectors 
Vi = h~ l diXi. Take \I/ G "Koo suited to (Xj(0); Vi(0)). Then by Proposition 
El there is a unique frame ip : Tj — ► !Koo adapted to the curve Xj with 
ip(0) = ^ this frame is defined as the unique solution of the differential 
equation 8(1]) = —Sie^ip with 5j given by eq. According to the proof of 

Proposition El the latter equation is equivalent to the system of equations 

diV k = PkA, Pki = ~(v k , dA), k^i. 

So, in order to determine the rotation coefficients (3ki : Tj 1— > M for 1 < k < N, 
k i, one has to solve the latter system of ordinary differential equations 
with the initial data 0^(0) = Ap(efc). Thus, we produced the functions hi 
and /3ki, or, what is equivalent, the Clifford elements 

>Sj 2 ^ ^ Pki^k hi^oc 
kjti 

for a given curve X, : CPj 1— > X. We say that hi, /3 ki are read off the curve Xj. 

Definition 11 The canonical discretization of the curve Xj : 3^ 1— ► % with 
respect to the initial frame \I/ G "Koo suited to (Xj(0); "Oj(O)) is the function 
X\ = A^(e ) : 1 1— > DC, where ip e : CP| 1— ► CKqo zs i/je solution of the discrete 
moving frame equation = — Ejej-0 6 




and t/ie initial condition if) e (0) = ^ . 

In other words, for the canonical discretization the data hi and j3ki are read 
off the continuous curve. Obviously, tp e is an adapted frame for the discrete 
curve X-. As e — > 0, the canonical discretization X\ converges to Xj with 
the rate 0(e) in C°°. 

Proposition 6 (approximation for C— surfaces). Let m = M = 2, and 

let there be given: 
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• two smooth curves Xi : Tj i— > ~R N (i = 1,2), intersecting orthogonally 
atX = X 1 (0) = X 2 (0), 

• a smooth function F : 3> 12 i— > K.. 

Assume \I/ G !Koo suited for (X; {^(O), ^(O)). Then, for some r > 0: 

1. There exists a unique C-surface x : -£>o( r ) l— * that coincides with 
Xi on Ti fl B (r) and satisfies 

difa - d 2 /3 21 = 2T. 

2. Consider the family of discrete C-surfaces {x e : Bq{t) — > M N } 0<t<tl 
defined as the solutions to the Goursat problem^ with the data 

v^ = vd, Ht = H { \^ BU = B ki \ n , r = r^ 2 , 

where the functions Hi and B ki are read off the smooth curves X i; so 
that x e \yc are the canonical e- discretizations of Xi. This family of 
discrete C-surfaces 0(e) -converges in C°° to x. 

Proof follows from Lemma 0] and Theorem EJ ■ 




Figure 6: Two curves determine a discrete C-surface via Proposition |3 



5.2 Goursat problem and approximation for an 
orthogonal system 

Theorem 6 (approximation of an orthogonal system). Let mim — 

l)/2 smooth C-surfaces §>ij : 7^ i— > M. N be given, labelled by 1 < i ^ j < 
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m with §>ij = Sji. Assume that for a given 1 < i < m all surfaces Sjj 

intersect along the curvature lines X, = S^- Set hi = [d^X^, assume 
that hi(0) 7^ 0, and set further Vi = K[ l diXi. All curves Xi intersect at one 

point X = Xl(0) = . . . = X m (0) orthogonally. Let ^ e "K^N) be suited 

to (X, Vi(0), ■ ■ ■ , v m (0)). Construct the discrete C-surfaces S^- : IPf^ i— »■ M N 

according to Proposition® with the functions = {difiij — dj(3ji)/2 for 
1 < i < j < m. Then for some r > 0: 

1. There exists a unique orthogonal system x : B (r) i— > ~R N , coinciding 
with §>ij on n B {r). 

2. There exists a unique family of discrete orthogonal systems {x e : Bq(t) h- > 
M Ar }o<e<e coinciding with Sf^ on CPf^- H Bq{t). The family x e converges 
to x with the rate 0(e) in C°° . 

Proof. C-surfaces are special two-dimensional conjugate nets, therefore the 
surfaces S^- can be supplied with the respective coefficients CV,-, Cji : 7^ i— > ffiL 
Now Theorem El can be applied. It yields the existence of a unique conjugate 
net x which has Xi as image of the z-th coordinate axes and Cij^Cji as 
coefficients on the respective 3V,. 

Similarly, discrete C-surfaces are special discrete two-dimensional con- 
jugate nets, so for each S^, the coefficients Cf^C^ are defined, and they 
converge in C°° to the coefficients (7^,(7^ of the respective S^-. Note that 
for any 1 < i < m the discrete surfaces Sf^ intersect along the discrete curve 
X- = §lj\j>i, which is the canonical discretization of the curve Xi. The con- 
clusion 2 of Theorem El is still valid if the Goursat data for cf^ in fK->2l) are 
taken as Cfj instead of CV,-, i.e. read off rather than S^. Thus, Theorem 
El delivers a family of discrete conjugate nets {x e } that is 0(e)-convergent 
to x. Since the restrictions of x e to the coordinate planes CPf - are discrete 
C-surfaces, the discrete variant of Lemma El implies that the nets x e are ac- 
tually discrete orthogonal systems. It remains to show that the limiting net 
x is with necessity an orthogonal one. But this follows immediately from the 
fact of C°°-convergence and eq. (f75j): indeed, we can conclude that for the 
net x everywhere holds 

Vjdi + Vii'j = <^> (di, Vj) = 0, 

that is, diX ■ djX = 0. ■ 

This result immediately implies Theorem^from Introduction: to bring its 
assumptions into the form of Theorem H above, one needs only to introduce 
a curvature line parametrization on surfaces 3^ such that the curves Xi are 
parametrized by arc-length. 
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5.3 Goursat problems and approximation for 
Ribaucour transformations 

Now we consider again a two-dimensional (M = 2) discrete orthogonal sys- 
tem described by eqs. (|1U5|) - ([1U8|) . but perform a different continuous limit, 
leading to two curves enveloping a circle congruence (Fig. HJ). In other words, 




Figure 7: A pair of curves enveloping a circle congruence 

we set m = m' = 1, ei = e, e2 = 1. Recall that in this case 

B £ (r) = B e (r) x {0, 1}, where B e (r) = T x n [0,r] = (eZ) n [0,r], 

so that only two lines in are considered. Therefore, a function on 7\ 2 
is conveniently considered as a pair of functions on 7\, or even as just one 
function on if one is interested in -Bg( r ) x {0} only (this will be the case 
for the function a below). We denote the shift r 2 by the superscript "+". 
Set 

021 = e« , P12 = W12 + e(N 2 p 21 - - a); (116) 

with a suitable function a : CPf i— > R. 

As the Banach space where the system lives we choose, exactly as before, 

X = , k{i)} x RW-^iP} x R 2 {h u h 2 } x R{a} . 

In the present case we have: 

iV2 = i + 0(e 2 ), iV 2 2 = l-i^/3 fc 2 2 , n =l- e -a(3 12 + 0{e 2 ), (117) 

where the constants in O-symbols are uniform on compact subsets of X. 
However, the system itself is now well defined not on all of X but rather on 
the subset where iVf > 0: 

© = fapMM*) : <4 } CI ( 118 ) 

k+2 
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Goursat Problem 4 (for a pair of curves enveloping a circle con- 
gruence). Given G functions HI : 7\ i— > K and -Bj^ : i— > R /or 
2 < k < N , and A e : 3^ i— > R ; as we// as iae rea/ numbers H\ and B e k2 for 
1 < k < N, k 2, /md a solution to the equations hlO/fy - HUIfy with hi Iffy 
on B e (r) with the Goursat data 

#)) = * e , M^=#i e , /W^k, 0^=^, 
h 2 (0) = HI, 13*2(0)= By. 

Lemma 5 Tae hyperbolic system hlO/fy - HUIfy with hi Iffy is consistent. 
Goursat problem^ satisfies condition (F) on the subset D. The limiting 
system for e = consists of equations \6Jj\) with hftffy for i — 1, h8(fy with 
(EIP fori = 2, and 





d\h 2 


= p 12 (h 1 + ah 2 /2) , 




(119) 




di(3 k2 


= Mf3 k i + af3 k2 /2) 


, k>2, 


(120) 






= 2(N 2 (3 21 - 6 - a) 


+ afi 2 /2, 


(121) 


hi 


-h 


= ah 2 , 




(122) 


til 


-Phi 


= a(3 k2 , k> 2, 




(123) 


til 


-hi 


= 2(N 2 a-p 21 ), 




(124) 



and describes a pair of continuous curves enveloping a circle congruence. 

Proof. Consistency is shown exactly as in Lemma HJ the limiting system is 
calculated directly from (jl()f)j) - (jl()8j) using ()117j) . One sees that the system 
of the lemma coincides with (|9"%|)-(fT0"3"j) for M = 2. * 

Proposition 7 (approximation of a pair of curves enveloping a circle 
congruence). Let m = 1, M = 2, and let there be given: 

• a smooth curve X : 7i t— > M. N , 

• a smooth function A : i— ► R ; 

• a point X + (0) G R N . 

Set H 2 (0) = \X+(0) -X(0\\, v 2 (0) = iJ 2 (0)- 1 (X+(0) - X(0)). Assume 
that \I/ G "Koo is suited for (X^O); #i(0)) ; and set Vk(0) = A$,(e k ) for k ^ 2. 
Then, for some r > 0: 

1. There exists a unique curve X + : ?i n B (r) i — ► R N through the point 
X + (0) such that the pair (X, X + ) envelopes a circle congruence, and 

\d 1 X + \-\d 1 X\= A-\X + - X\. (125) 
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2. Consider the family of pairs of discrete curves (X e , (X e ) + ) 0<e<£1 defined 
as the solutions to the Goursat problem^ with the data 

= H\ = Hi \ye , B% x = B kl \ n , A e — A\je, 
HI = H 2 (0), B e ka = -2(v 2 {0),v k (0))- 

where the functions Hi and B k i are read off the smooth curve X, so 
that X e is the canonical e- discretization of X. These discrete curves 
0(e) -converge in C°° to (X, X + ). 

Proof follows from Lemma El and Theorem [2j ■ 

Theorem 7 (approximation of a Ribaucour pair). Let, in addition 
to the data of Theorem there be given m curves X^ : O 5 , i— ► M. N with 
a common intersection point X + = Xi (0) = . . . = X+(0) ; and such that 
each pair (Xi,X^) envelopes a circle congruence. In addition to the discrete 
surfaces S|- from Theorem® construct discrete curves (Xl) + according to 
Proposition^ with the functions Ai = (\diX^\ — \diXi\) /\Xf — X{\. Then 
for some r > 0: 

1. There exists a unique Ribaucour pair of orthogonal systems x,x + : 
Bo{r) \— > M> N such that x coincides with §>ij on 7^ n Bq(t), and x + 
coincides with Xf on CPj fl B (r). 

2. There exists a unique family of Ribaucour pairs of discrete orthogonal 
systems {x e : B e (r) \— > R Ar } < e < eo coinciding with S^- on (O 5 ^- x {0}) fl 
B(r) and coinciding with (Xf) + on (7\ x {1}) fl B(r). The family x e 
converges to the pairs x,x + with the rate 0(e) in C°° . 

Proof is similar to that of Theorem® with the only change in the conculding 
argument: for the limiting Jonas pair x, x + of orthogonal nets we derive from 
eq. (JZSD: 

VMvf + ViV M =0 =$> vf, V M ) = 0, 

which is the defining property of the Ribaucour pair. ■ 

Letting not one but two or three directions of the orthogonal system 
remain discrete in the continuous limit (so that w! = 2 or ml = 3), one 
arrives at the following statement. 

Theorem 8 (permutability of Ribaucour transformations). 
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1. Given an m- dimensional orthogonal system x(-,0, 0) : B (r) i— > K.^ 
and zts two Ribaucour transformations x(-,l,0) : -£>o(r) i— > and 
x(-,0, 1) : i?o( r ) l— * ^ N , there exists a one-parameter family of orthog- 
onal systems x(-, 1, 1) : Bo(r) i— > K ina£ are Ribaucour transforma- 
tions of both x(-, 1, 0) and 0, 1). Corresponding points of the four 
conjugate nets are concircular. 

2. Given three Ribaucour transformations 

x(; 1, 0, 0), X(; 0, 1, 0),X(; 0, 0, 1) : B (r) i-> R N 

of a given m-dimensional orthogonal system x(-, 0, 0, 0) : Bo(r) i— > M^, 
as well as three further orthogonal systems 

X(; 1, 1, 0), X(; 0, 1, 1), X(; 1, 0, 1) : B (r) h-> 

swca taai 1, 1, 0) is a Ribaucour transformation of both x(-, 1, 0, 0) 
and #(-, 0, 1, 0) etc., taere exists generically a unique orthogonal system 

x{-,l,l,l): B {r)^M N 

which is a Ribaucour transformation of all three 1, 1, 0), 0, 1, 1) 
and x(-, 1, 0, 1). 



5.4 Example: elliptic coordinates 

The simplest nontrivial example, to which the above theory can be applied, is 
the approximation of two-dimensional conformal maps by circular patterns. 
Starting with a conformal map F : M 2 — > ffi 2 , i.e. M = N = 2, 

h = \diF\ = \d 2 F\ and d x F ■ d 2 F = 0, 

one calculates the metric and rotation coefficients according to equations (}4"3)) 
and (|3HJ) . and the quantity 7 according to pi4jl : 

hi = h 2 = h, (3\2 = ^— , P21 = -r 1 , 7 = di[3 12 = -d 2 (3 21 

n h 

To construct a discrete approximation of F, solve the Goursat problem |3] 
for e > with the data B\, and T £ that are close to the values of the 
respective functions hi, flij and 7 at the corresponding lattice sites. A good 
choice is, for example, to prescribe 

= + e/2,0), = h(0,£ 2 + e/2), 

SJi(£i) = + e/2,0), £f 2 (&) = ^(0,6 + e/2), 

r e (£i,6)=7(a + e/2,6 + e/2), 
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Figure 8: Approximation of elliptic coordinates, e = n/10 and e = n/20. 



for £1, £ 2 = i€, < i < R. If ^ e : B t (eR) — > IKoo is the frame of the respective 
solution to the system (jl05|) - (jl08|) . then the function 

F e : B\Re) h- R 2 , F e (£ 1; £ 2 ) = 7r(^(£i, ^"WCCi, £2)) 

is a two-dimensional discrete orthogonal system, i.e., the points -^(^i,^), 
+ e, £2), -^ e (£i,^2 + e) and -F e (£i + e, £2 + e) lie on a common circle 
C e (£i,£ 2 ) in M 2 , and F e (&,6) = + 0(e). 

This is illustrated with the planar elliptic coordinate system: 

^(£i,£ 2 ) = (cosh(6)cos(6),sinh(ei)sin(6)), 
whose coordinate lines are ellipses and hyperbolas. One finds 

/ i =(sinh 2 (6) + sin 2 (6)) 1/2 , 
Pn = sinh(26)/(2/i 2 ), fa = sin(2£ 2 )/(2/i 2 ), 
7 = (l-cosh(2ei)cos(26))/(4/i 4 ). 
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Results are displayed in Fig. |HJ Their left sides show the original coordinate 
lines of F, and on the right sides the circles C e (£i,£2) are drawn; each inter- 
section point of two coordinate lines on the left corresponds to an intersection 
point of four circles on the right. There are small defects (the circles do not 
close up) on the very right of the pictures since, in contrast to F, the discrete 
maps F e are not periodic with respect to £ 2 - 
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